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(57)Abstract: 

PROBLEM TO BE SOLVED: To provide a 
fluid flow measuring technology allowing the 
accurate measurement of fluid flow in a 
complicated flow field from a remote place. 
SOLUTION: The system for measuring fluid 
flow comprises: a laser oscillator 1 1 ; a laser 



sheet forming scanning optical system 13 
for irradiating the fluid flow field 14 with laser beams oscillated from the oscillator 
1 1 , in a sheet shape; an image pickup means 29 for picking up a 
two-dimensional particle trace image on a laser sheet 17; a timing control means 
30 for taking timing to synchronously drive the laser oscillator 1 1 and the image 
pickup means 29; and an image processing means 35 for processing the 
picked-up particle trace image. The image processing means 35 compares and 
analyzes the luminance pattern distribution of the particle images at two different 
times to measure the moving direction and distance of individual particles, thus 
measuring fluid flow and flow velocity distribution in the fluid flow field. 
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CLAIMS 
[Claim(s)] 

[Claim 1] The laser oscillation equipment which oscillates a laser beam, and the 
scan optical system for laser sheet formation which makes the oscillated laser 
beam supply to the floating in the hall of a fluid in the shape of a sheet, An image 
image pick-up means to picturize the two-dimensional particle track image on 
the laser sheet from this scan optical system, A timing control means to take and 
carry out synchronization actuation of said laser oscillation equipment and the 
image image pick-up means for timing, The brightness pattern of the particle 
track image of 2 time of day picturized with said image image pick-up means is 
compared and analyzed. The floating instrumentation system of the fluid 
characterized by constituting so that it may have an image-processing means to 



measure each migration direction and movement magnitude of a particle and 
floating and the velocity distribution of a fluid may be measured from the 
movement magnitude of each particle by this image-processing means. 
[Claim 2] Said image image pick-up means is the floating instrumentation system 
of the fluid [ equipped with the picture transmission means which carries out 
optical transmission of the two-dimensional particle track image on the laser 
sheet from scan optical system, and an image pick-up means to picturize the 
two-dimensional particle track image by which optical transmission was carried 
out ] according to claim 1. 

[Claim 3] Said a majority of picture transmission means bundle the optical fiber 
of a book, and are made to unify, and it consists of flexible image guides to 
which flattening of the ends side was carried out. Said image pick-up means It is 
the floating instrumentation system of the fluid according to claim 2 constituted 
by consisting of CCD cameras so that said image guide might disassemble into 
each pixel the image by which image formation was carried out to one fiber end 
face with the objective lens with an optical fiber and the same image might be 
transmitted to the other end side by the side of a CCD camera. 
[Claim 4] said scan optical system and an image image pick-up means are 



combined in one, and are constituted — having — the floating field of a fluid 
— dimensional [ 1 ] or two-dimensional — the floating instrumentation system of 
the fluid according to claim 1 constituted movable. 

[Claim 5] Said scan optical system is the floating instrumentation system of the 
fluid according to claim 1 which established the picture transmission means 
which forms in the floating in the hall of a fluid the laser sheet which intersects 
perpendicularly, and carries out optical transmission of the two-dimensional 
particle track image on this laser sheet that intersects perpendicularly, 
respectively. 

[Claim 6] Irradiate the laser beam oscillated from laser oscillation equipment in 
the shape of a sheet at the floating field of a fluid, and a laser sheet is formed. 
While picturizing the particle track image on this laser sheet with an image image 
pick-up means and comparing and analyzing the brightness pattern of the 
particle track image of 2 picturized time of day Carry out averaging processing of 
the particle track image of the picturized time series, and it asks for the 
brightness pattern of an optical fiber array image. The floating measurement 
approach of the fluid characterized by calculating each migration direction and 
movement magnitude of a particle by removing the brightness pattern of this 



optical fiber array image from the brightness pattern of said particle track image, 
and measuring floating of the fluid in said fluid flow field. 

DETAILED DESCRIPTION 

[Detailed Description of the Invention] 
[0001] 

[Field of the Invention] This invention relates to the floating instrumentation 
system and its measurement approach of the fluid which measures the rate of 
flow and flow direction of a fluid which start the floating measurement technique 
of the fluid which measures floating of the fluid in a complicated floating field to 
high degree of accuracy and a precision, especially flow the inside of a closed 
space. 
[0002] 

[Description of the Prior Art] It is very important for high-degree-of-accuracy-izing 
of the plant behavior prediction in a nuclear power plant, a thermal power station 
plant, a chemical processing plant, etc., a plant diagnosis, and a plant 
performance evaluation and device amelioration to measure floating of the fluid 
in a complicated floating field to high degree of accuracy and a precision. The 



floating measurement of the floating measurement technique of the fluid under 
the severe environment which differs in the plant especially the external world, 
and the environment of system level is in a very difficult situation because of fluid 
flow measurement conditions with complicated flow. 

[0003] Especially, the fluid flow measurement of the rate of flow of a fluid, a flow 
direction, etc. in the heat floating field under a severe environment is very 
difficult including access nature because of complicated fluid flow rate 
measurement conditions like [ in a reactor pressure vessel or the heat 
exchanger of a thermal power station plant ]. 

[0004] However, research of the visualization technique in a heat floating field 
progresses in recent years, and a particle image current meter (it is called PIV 
below Particle Image Velocimetry:.) measurable to high degree of accuracy and 
a precision is being developed in floating of the fluid in a complicated floating 
field. By using together the image guide which bundled this PIV and optical fiber, 
floating measurement of the fluid in the external world and the closed space 
which differs in an environment is expectable like [ in a reactor pressure vessel ]. 
[0005] The fiberscope which transmits an image directly on the other hand using 
the image guide which bundled the optical fiber is developed as a 



medical-application endoscope. The image guide was rich in flexibility\and has 
achieved the duty important as a picture transmission means for an observation 
means. 

[0006] The image guide developed by medical application is applied also to the 
industrial field, and is applied to the interior check of an aircraft engine, and 
observation of a short pipe or a reactor reactor core by the end of today. 
[0007] Moreover, by technological innovation in recent years, the conventionally 
difficult long-distance picture transmission also becomes possible by 
development of the small optical fiber of optical loss, and it is in the limelight as a 
new maintenance inspection means from the remote place in a large-scale plant 
facility. 

[0008] On the other hand, the research on PIV has improvement in time 
resolution and spatial resolving power, reservation of fluid flow measurement of 
high degree of accuracy, three-dimensions measurement, etc. in the trend 
performed briskly from both sides of hardware and software. One of the 
researches-and-developments trends of these PIV is the movement toward 
"improvement in precision" of the fluid rate-of-flow data obtained by PIV, and 
other one is the movement toward "applicability amplification of PIV." 



[0009] 

[Problem(s) to be Solved by the Invention] Although PIV applied to the floating 
measurement technique of a fluid has been developed as a laboratory mold, it is 
requested that the object for development of this PIV technique should be 
expanded from a laboratory to utilization. The development request to utilization 
of the fluid flow measurement to real fluidization level is strong also in the ED of 
a direction to which the floating measurable range of the fluid by the PIV 
technique is expanded. 

[0010] This invention was made in consideration of the situation mentioned 
above, and aims at offering the floating instrumentation system and its 
measurement approach of the fluid which can measure floating of the fluid in a 
complicated floating field with a sufficient precision to accuracy. 
[001 1] Other objects of this invention are to offer the floating instrumentation 
system and its measurement approach of the fluid which can measure the rate 
of flow and flow direction of a fluid of a floating field in the external world and the 
closed space which differs in an environment with a sufficient precision from 
remoteness to accuracy. 

[0012] The object of further others of this invention optimizes the algorithm which 



performs an image processing by the particle track pursuing method (Particle 
streak tracking velocimetry) for having used the image guide, and detects the 
floating condition of a fluid, and is to offer the floating instrumentation system 
and its measurement approach of the fluid which can perform floating 
distribution of a fluid with a sufficient precision to accuracy. 
[0013] 

[Means for Solving the Problem] In order that the floating instrumentation system 
of the.fluid concerning this invention may solve the technical problem mentioned 
above The laser oscillation equipment which oscillates a laser beam as indicated 
to claim 1 , The scan optical system for laser sheet formation which makes the 
oscillated laser beam supply to the floating in the hall of a fluid in the shape of a 
sheet, An image image pick-up means to picturize the two-dimensional particle 
track image on the laser sheet from this scan optical system, A timing control 
means to take and carry out synchronization actuation of said laser oscillation 
equipment and the image image pick-up means for timing, The brightness 
pattern of the particle track image of 2 time of day picturized with said image 
image pick-up means is compared and analyzed, and it has an 
image-processing means to measure each migration direction and movement 



magnitude of a particle, and it constitutes so that floating and the velocity 
distribution of a fluid may be measured from the movement magnitude of each 
particle by this image-processing means. 

[0014] Moreover, in order to solve the technical problem mentioned above, it 
sets to the floating instrumentation system of the fluid of this invention. As 
indicated to claim 2, said image image pick-up means As it has the picture 
transmission means which carries out optical transmission of the 
two-dimensional particle track image on the laser sheet from scan optical system, 
and an image pick-up means to picturize the two-dimensional particle track 
image by which optical transmission was carried out and being further indicated 
to claim 3 Said a majority of picture transmission means bundle the optical fiber 
of a book, and are made to unify, and it consists of flexible image guides to 
which flattening of the ends side was carried out. Said image pick-up means It 
consists of CCD cameras, and said image guide is constituted so that an optical 
fiber may decompose into each pixel the image by which image formation was 
carried out to one fiber end face with the objective lens and the same image may 
be transmitted to the other end side by the side of a CCD camera. 
[0015] Furthermore, in order to solve the technical problem mentioned above, it 



sets to the floating instrumentation system of the fluid of this invention. As 
indicated to claim 4, said scan optical system and an image image pick-up 
means it combines in one and constitutes -- having - the floating field of a fluid 
—dimensional [ 1 ] or two-dimensional - said scan optical system, as it was 
constituted movable and indicated to claim 5 The laser sheet which intersects 
perpendicularly is formed in the floating field of a fluid, and the picture 
transmission means which carries out optical transmission of the 
two-dimensional particle track image on this laser sheet that intersects 
perpendicularly is established, respectively. 

[0016] In order that the floating measurement approach of the fluid concerning 
this invention may solve the technical problem mentioned above As indicated to 
claim 6, irradiate the laser beam oscillated from laser oscillation equipment in the 
shape of a sheet at the floating field of a fluid, and a laser sheet is formed. While 
picturizing the particle track image on this laser sheet with an image image 
pick-up means and comparing and analyzing the brightness pattern of the 
particle track image of 2 picturized time of day Carry out averaging processing of 
the particle track image of the picturized time series, and brightness pattern 
distribution of an optical fiber array image is searched for. It is the approach of 



calculating each migration direction and movement magnitude of a particle, and 
measuring floating of the fluid in said fluid flow field by removing the brightness 
pattern of this optical fiber array image from the brightness pattern of said 
particle track image. 
[0017] 

[Embodiment of the Invention] One operation gestalt of the floating 
instrumentation system of the fluid concerning this invention and its 
measurement approach is explained with reference to an accompanying 
drawing. 

[0018] Drawing 1 is principle drawing showing the configuration of the floating 
instrumentation system of the fluid concerning this invention. 
[0019] The floating instrumentation system 10 of the fluid of this invention is 
equipped with the laser oscillation equipment 1 1 which oscillates the laser beam 
of high power as the light source. To laser oscillation equipment 11, it is small, 
and semiconductor laser, a pulse laser, etc. which oscillate a high power laser 
beam are used. For example, semiconductor laser with a main oscillation 
wavelength of 808nm is used as the light source with the continuation light of 
output 16W. 



[0020] The laser beam oscillated from laser oscillation equipment 11 is guided 
through the fiber 12 for light transmission at the scan optical system 13 for laser 
sheet formation. The scan optical system 13 makes the laser beam oscillated 
from laser oscillation equipment 1 1 project on the floating field 14 of a fluid in the 
shape of a sheet. 

[0021] The floating field 14 of a fluid is formed in the external world and the 
closed space which differs in an environment. Moreover, according to the case, 
the floating field 14 of a fluid may be formed with the gestalt opened in the 
external world. Inside [ the floating fields 14 are a gas and a liquid ] consists of 
one of fluids at least. It is aimed concrete at the rate field of the fluid in the closed 
space in the heat exchanger of the downcomer section of a reactor pressure 
vessel, a reactor core shroud, and a thermal power station plant, or a steam 
generator etc. Moreover, the floating field of the fluid [ structured division / in 
contact with the stream or air current of a marine vessel, the aircraft, or an 
automobile ] of the circumference of equipment and materials can also be made 
into an object. In addition, principle drawing shown in drawing 1 shows the 
example of an experiment of the floating field of water as a floating field 14 of a 
fluid. In a container 15 or piping, from the nozzle 16, water is gushed and the 



floating field 14 of the water as a fluid is formed. 

[0022] The scan optical system 13 irradiates the oscillation laser beam from 
laser oscillation equipment 1 1 in the shape of a sheet at the floating field 14 of a 
fluid, and forms the laser sheet 17 in the fluid flow field 14. The floating field 14 of 
a fluid is visualized with this laser sheet 17. In the floating instrumentation 
system 10, in order to secure the spatial resolving power of the depth direction of 
the fluid flow field 14, the oscillation laser beam is irradiated in the shape of a 
sheet. 

[0023] The scan optical system 13 is the laser beam exposure optical system for 
creation of the laser sheet 17, as shown in drawing 2 . If semiconductor laser is 
used as laser oscillation equipment 1 1 , outgoing radiation of the oscillation laser 
beam from semiconductor laser will be carried out with a big angle of divergence 
from the fiber 12 for light transmission. The laser beam by which outgoing 
radiation was carried out does not serve as light of axial symmetry, but has an 
astigmatism discrepancy, and configuration conversion of a laser beam 
becomes complicated. 

[0024] For this reason, the lens group 18 as an optical means for laser sheet 
formation for forming the laser sheet 17 from the laser beam by which outgoing 



radiation was carried out from the fiber 12 for light transmission is used for the 
scan optical system 13. The optical means for laser sheet formation is 
constituted combining the collimate lenses 19a and 19b of two sheets with which 
this lens group 18 makes a concurrency laser beam collimate an outgoing 
radiation laser beam, the long focal lens 20 which carries out the waist of the 
laser beam which became concurrency light with collimate lenses 19a and 19b 
to about about 1-2mm in the thickness direction of the laser sheet 17, and the 
cylindrical lens 21 which forms the laser sheet 17. A plano-convex lens is used 
for collimate lenses 19a and 19b. 

[0025] The laser sheet 17 is formed in the floating field 14 of a fluid of the laser 
beam irradiated in the shape of a sheet from the scan optical system 13, and the 
floating field 14 of a fluid is visualized in the shape of a sheet. And 
measuring-range 17a of the laser sheet 17 is countered, and the picture 
transmission means 22 is established. As shown in drawing 1 , the windowpane 
23 for transparency is installed in container 15 wall near the picture transmission 
means 22 so that the image pick-up of the image on the laser sheet 17 may be 
attained with the image image pick-up means 29. 

[0026] The image image pick-up means 29 has CCD camera 24 as a picture 



transmission means 22 to transmit the image on the laser sheet 17, and an 
image pick-up means to picturize a two-dimensional particle track image. The 
picture transmission means 22 makes the two-dimensional particle track image 
on the laser sheet 17 transmit to CCD camera 24 as an image pick-up means. 
As the picture transmission means 22 is shown in drawing 3 , the flexible image 
guide 26 which bundled thousands of to tens of thousands of optical fibers 25, 
and was made to unify (melting) is used. As shown in drawing 3 , it is installed 
and an image guide 26 is bundled so that each location of the ends side of the 
bundled optical fiber 25 may correspond to accuracy. A flat surface is made to 
the ends side of an image guide 26, it disassembles into the optical fiber 25 as 
each pixel fiber the image by which image formation was carried out to one fiber 
end face with the objective lens 27, and is transmitting the same image to the 
fiber end face by the side of CCD camera 24. It is recorded on a CCD camera 
through the camera lens 28 from a fiber end face. The image image pick-up 
means 29 is made to carry out the direct image pick-up of the two-dimensional 
particle track image on the laser sheet 17, without having an image guide. 
[0027] How an image appears by deciding the limitation (resolution) of the image 
of an image guide 24 by the diameter of the single optical fiber 25 and the 



method of an array is decided by the total number and arrangement of the 
optical fiber 25 bundled by the image guide 26. Optical distribution of the 
quantity of light smaller than the core diameter of an optical fiber 25 is equalized 
while an optical fiber 25 is transmitted, and in each core cross section, 
brightness becomes uniform. Therefore, as for an image guide 26, the core 
diameter of an optical fiber 25 serves as a smallest unit of the spatial resolving 
power as a pixel (pixel). 

[0028] As for the actual image guide 26, many optical fibers 25 of a book are 
arranged tidily and densely with ****** (roppo dense array) structure. If the 
grid-like body of spacing almost equal to optical fiber spacing is observed by this 
image guide 26, a moire interference fringe unrelated to the intensity distribution 
of an objective image will arise. 

[0029] Moreover, in the image guide 26 by which welding of the clad part is 
carried out between optical fibers 25, if the thickness of a clad becomes thin, the 
case where the light included in the core of one optical fiber 25 leaks to the next 
core will arise. For this reason, the boundary of an image spreads and the 
contrast of an image falls. For this reason, it becomes conditions that an image 
guide 26 does not have the too thin thickness of a clad when roppo dense array 



structure is adopted. 

[0030] On the other hand, the floating instrumentation system 10 of a fluid is 
equipped with the timing control means 30. The timing control means 30 
consists of a timing scheduler 31 and a synchronizer 32. While laser oscillation 
equipment 11 is made to oscillate to the timing from the timing scheduler 31, a 
synchronization is taken by this oscillation timing and synchronizer 32, and CCD 
camera 24 is made to drive. That is, the timing control means 30 takes the 
synchronization with laser oscillation equipment 11 and CCD camera 22 as an 
actuation means, and is made to drive it. 

[0031] A thing (resolution [ which has sensibility in a near infrared ray field ] VGA 
(640x480pixels), and frame rate 30Hz, and 8 bits (monochrome 256 gradation)) 
is used for a CCD camera. The image processing of the picture signal of the 
analog which CCD camera 24 took the oscillation timing from laser oscillation 
equipment 11 and a synchronization, was picturizing the two-dimensional 
particle track image on the laser sheet 17, and was picturized with CCD camera 
24 is sent and carried out to the image-processing means 35, and the velocity 
distribution and flow direction of a fluid in the floating field 14 of a fluid are 
measured by the particle track pursuing method mentioned later. The camera of 



a read-out method (progressive scan) is serially used for CCD camera 24, and 
the charge storage time of odd lines of arranged CCD and even lines is 
simultaneously performed to it, for example. This CCD camera 24 is excellent in 
vertical resolution compared with the usual camera. 

[0032] Moreover, a digital picture signal is sent through the signal lines 37, such 
as a PCI bus, from the frame grabber board 36 and this frame grabber board 36 
as an A/D converter which digitizes the analog picture signal from CCD camera 
24, and the image-processing means 35 has the computer 38 which carries out 
digital image processing. 

[0033] The floating instrumentation system 10 of this fluid carried out analog 
photography with CCD camera 24, and carried out digital processing with the 
image-processing means 35, and the analog photography and the 
digital-recording method which carries out digital recording are used for it. 
[0034] The floating instrumentation system 10 of the fluid shown in drawing 1 
showed the example corresponding to a laboratory. On the occasion of 
utilization, as shown in drawing 4 and drawing 5 , the scan optical system 13 for 
laser sheets and the image image pick-up means 29 are unified. The picture 
transmission means 22 of the image image pick-up means 29 united with the 



scan optical system 13 consists of drawing 4 as a probe 40 for visualization 
which can view laterally. The 45-degree rectangular prism 43 carries out picture 
transmission of the optical path by the fiberscope 44 which constituted the 
bending picture transmission means from the optical aperture (windowpane for 
transparency) 41 installed at probe 40 head. Drawing 4 (A) shows the example 
which used the floating instrumentation system of this unified fluid for the floating 
field of the water as a fluid in the container 15 of drawing 1 . In drawing 5 , the 
unified scan optical system 13 and the unified picture transmission means 22 of 
the image image pick-up means 29 are stored in the body casing 45. 
[0035] Next, the floating measurement approach of a fluid is explained. 
[0036] In the floating instrumentation system 10 of the fluid shown in drawing 1 , 
the comparatively long luminescence time amount of the semiconductor laser 
used for lighting is taken, and the pathline (locus) of a particle is picturized to 
CCD camera 24. And it considers that the particle track image of this 2 picturized 
time of day that is different minute time in time is a brightness pattern, the 
computer 38 of the image-processing means 35 performs the image processing 
of a comparison and analysis, and the movement magnitude of each particle is 
presumed. In the floating metering device 10 of this fluid, this technique is called 



the particle track pursuing method (Particle streak tracking velocimetry). This 
particle track pursuing method can give the information on enough particle 
images to the optical fiber single fiber used as the minimum spatial resolving 
power, without enlarging particle **** by recording the pathline of a particle 
beyond the need. 

[0037] The particle track image expresses the image of the pathline in the time of 
day which floats in the water on the laser sheet 17 (fluid), and has each particle 
(foreign matters, such as suspended matter) which carries out diffusion 
distribution, and is constituted by the pixel (pixel) with the digitized brightness. In 
process of the image processing by the image-processing means 35, rectangle 
field 17a (refer to drawing 1 ) limited among each particle track image is taken 
out, and an image processing is carried out. A certain value of one point in the 
digitized particle track image is defined as a "brightness value", and that from 
which this brightness value is distributed over the field of rectangle field 17a is 
called a "brightness pattern." 

[0038] The algorithm which presumes the movement magnitude of each particle 
in water (fluid) from a particle track image judges the obtained particle track 
image to be change of luminance distribution, and analyzes quantitatively the 



amount which moves into time amount with the brightness pattern of the particle 
on the laser sheet 17 about each particle. A cross-correlation method is used for 
the analysis of the movement magnitude of this particle. 
[0039] However, not only the brightness pattern of a particle but the brightness 
pattern of the array (bag structure to stack) of the optical fiber 25 of an image 
guide 26 exists in the particle track image picturized by CCD camera 24 through 
the image guide 26. The array brightness pattern of an optical fiber 25 serves as 
evil in the case of the analysis by the cross-correlation method, and causes an 
incorrect response (mistake vector) to the brightness pattern of a particle. 
[0040] Then, in the floating instrumentation system 10 of the fluid concerning this 
invention, in order to attain optimization of floating measurement of the fluid by 
the particle track pursuing method, array clearance processing (background 
processing) of the optical fiber 25 used as the hindrance of cross-correlation 
count is performed. 

[0041] the floating instrumentation system and its measurement approach of the 
fluid concerning an application-for-patent No. 207476 [ 2000 to ] official report - 
a particle image flow pattern measurement (PIV) ~ law - applying - this PIV - a 
freezing image is used in law. By the usual PIV method, if particle **** is smaller 



than one optical fiber 25 single fiber of an image guide 26, since the information 
on the already lost particle is not recoverable, the correlation coefficient 
computed is strongly influenced to a ************ pattern, serves as very serious 
evil at speed detection, and causes an incorrect response. 
[0042] moreover, PIV which will estimate the movement magnitude of a particle 
group if particle **** is made larger than one single fiber of the optical fiber 25 
which gives the minimum spatial resolving power — in law Since it is made 
desirable for there to be 6-8 particles in a reference aperture image also at the 
lowest, A comparatively large reference aperture image must be taken, and 
computational complexity increases enormously and will become not realistic as 
an algorithm (if the fields to search increase in number, computational 
complexity will increase in proportion to the 4th power). 
[0043] for this reason, PIV using a freezing image - in law, the parameter 
(particle ****, number density) demanded receives a limit dramatically. So, by the 
particle track pursuing method used by this invention, the comparatively long 
luminescence time amount of the semiconductor laser used for lighting was 
taken, the algorithm which estimates the movement magnitude of a particle to be 
each instead of record and a particle group for the pathline of a particle was built, 



and it mounted in software. The information on enough particle images can be 
given to the single fiber of the optical fiber 25 used as the minimum spatial 
resolving power, without taking large particle **** by recording the pathline of 
each particle beyond the need. The technique of this particle track pursuing 
method (Particle streak tracking velocimetry) becomes possible by using the 
semiconductor laser which can control the luminescence time amount of laser 
easily. 

[0044] By the particle track pursuing method used by this invention, calculation 
of the movement magnitude of a particle is mainly performed through six steps 
of the following procedures also including this background processing. 
[0045] (1) Clearance of a noise (background processing) 
The particle track image picturized with CCD camera 24 through the image 
guide 26 is processed, the addition average value of a brightness value is 
calculated for the particle track image obtained by time series for every pixel, 
and brightness pattern clearance of the fiber array of an optical fiber 25 is 
performed by lengthening an addition average value from the brightness value of 
each particle track image. 

[0046] The image-processing algorithm of brightness pattern clearance of an 



optical fiber array is explained below. 

[0047] When setting the particle track image of measurement data to x, the 
particle track image x is [Equation 1]. 

x= e +n ■ ( 1 ) 

It is come out and expressed. 

[0048] That is, in the particle track image x of measurement data, I think that the 
image with which the noise (brightness pattern of a particle) lapped with the 
signal (brightness pattern of an optical fiber array) is measured. Brightness 
pattern clearance processing of the optical fiber array by taking averaging is 
effective technique, when an SN ratio is small and it is possible to repeat 
measurement repeatedly on the case where there is no great difference also in 
the configuration frequency of a signal and a noise, or the same conditions. 
[0049] When the load of the time series element is carried out and a formula 1 is 
displayed, it is expressed like a formula 2. 
[0050] 
[Equation 2] 



x 1 (k) = s 1 (k) +n* (k) (2) 

i = i#a<Z>*ft^®#, 

k : ife@©fi^®ffcF*J©k#@(B t?£ -fe/Vlfgil. 

It is each 

[0051] In photography of a particle track image, the averaging x to the particle 
track image of M sheets, then M sheets (k) can express the number of sheets of 
the captured particle track image as follows. 
[0052] 
[Equation 3] 

[0053] The brightness pattern s of an optical fiber array is still the brightness 
value of a basis, when averaging (it ** by M) processing which becomes one M 
times the brightness value of this by M addition, and averages this is performed, 
since the same pattern will appear in a continuous particle track image, if it is the 
same experiment conditions (a light transmission system and light-receiving 
system). That is, the image in which only the optical fiber array in which a particle 
does not exist was reflected by averaging processing can be obtained, and if the 



brightness value (brightness pattern) of an optical fiber array is lengthened for 
every pixel from the particle track image with which an optical fiber array exists, it 
will become removable [ the brightness pattern of an optical fiber array ]. 
[0054] Since the background luminance of the acquired image becomes 
homogeneity mostly over the whole screen by this background processing, also 
not only in validity but in binarization processing, it becomes an effective means 
to correlation computation. 

[0055] (2) Perform binarization and an expansion-ized image processing as 
pretreatment for discriminating a particle from binarization and an 
expansion-ized image-processing particle track image. The locus image of the 
particle which floats in water (fluid) is influenced of the clad of an optical fiber 25, 
and will not become continuous, but will be recorded as an intermittent image 
(refer to drawing 6 (a)). Even if it carries out binarization processing of this 
intermittent particle image, the part in which information was lost by the clad is 
not recovered (refer to drawing 6 (b)). In order to make the locus image of these 
divided particles recognize as a migration locus of the same particle, an 
expansion-ized image processing is performed. An expansion-ized image 
processing will change all the pixels of the perimeter into 255 from 0, if the pixel 



of attention exists (expansion processing of 8 connection). With this algorithm, 
the intermittent particle image divided by repeating this expansion-ized 
processing twice is recovered as a continuous locus image (refer to drawing 6 
<c»- 

[0056] (3) Since the particle image of plurality (large number) exists in labeling 
and the particle track image of one detection of a particle image location, 
perform labeling which makes each recognize as each particle. Another number 
is given to a connected component which is different in the same number in each 
connected pixel. It calls it labeling to assign a number to such each field. A 
center position is computed by the particle image identified by labeling returning 
to a former image. A main coordinate is searched for by the following methods of 
elastic center (Particle image centroid). 

[0057] When the center of gravity which calculates X0 and Y0, and i and j are 

made into the coordinate of each point and I is made into the brightness value of 

the coordinate, centers of gravity X0 and Y0 are [Equation 4]. 
x m lLIM. {4 ) 



It is come out and expressed. The same is said of Y0. 



[0058] (4) Since the center position of a particle track image was detected by 
particle image pretreatment of calculation [ of the integer movement magnitude 
by correlation count ] (1) - (3), calculate the movement magnitude of each 
particle with a cross-correlation method. 

[0059] The particle track image which carried out digital processing of the 
two-dimensional particle track image in the 1st time of day picturized with CCD 
camera 24 of the floating instrumentation system 10 of a fluid In addition, a 
"reference image", The rectangle field limited among the "retrieval image" and 
the "retrieval image" in the particle track image which carried out digital 
processing of the two-dimensional particle track image in the 2nd time of day 
when a "reference aperture image" differs the rectangle field limited among 
"reference images" minute time from the 1st time of day is defined as a "retrieval 
aperture image", respectively. 

[0060] Movement magnitude presumption of the particle by the cross-correlation 
method is the technique of having the correlation value R expressed with a 
degree type being known between the "retrieval aperture image" and the 
"reference aperture image", and asking for the similarity of the brightness pattern 
of each particle, carrying out comparison examination with this correlation value 



R, and calculating the movement magnitude of each particle between 2 images. 

[0061] 

[Equation 5] 
R(X',Y',S,r]) 



J j(W-+<.y'+i)-/r)*-(/ a (jr'+i+f.y' + ;+^-/r) 



*-l mi-l 



(5) 



[0062] 
[Equation 6] 



/«v%r 0 (T rove O 0 

1 I > 2 = 



rt-1 m-1 

(6) 



[0063] As shown in a formula 5 and a formula 6, the correlation value R takes the 
values from -1 to 1 by what (normalization) is **(ed) with the mean square of the 
brightness value of a retrieval aperture image and a reference aperture image. 
Two images mean that the similarity between aperture images is so large that 
the value of the correlation value R is large completely in accordance with the 
time of the correlation value R being 1. 

[0064] In a formula 5 and a formula 6, 11 and 12 express the brightness of each 



pixel of a reference aperture image and a retrieval aperture image, and xi and 
eta are expressed in the relative location of a retrieval aperture image and a 
reference aperture image. The locations xi and eta where a correlation value 
serves as max are equivalent to movement magnitude deltaX' of the particle 
within an image, and delta Y 1 . Magnitude nxm of a reference aperture image and 
the relative location with a reference aperture image are determined from the 
min and maximum velocity which are predicted. 

[0065] (5) Since the integral value of the movement magnitude of a particle 
image can be found by correlation count of calculation (4) of subpixel movement 
magnitude, estimate subpixel movement magnitude by normal-distribution 
approximation (Gaussian peak-fit) used near [ three ] the peak of a correlation 
coefficient. It is independently computed about x and the direction of y, 
respectively. 

[0066] (6) Since the movement magnitude of a particle track image can be found 
per pixel (pixel) with the procedure more than a calibration, it is changed into the 
movement magnitude in real space by calibration. 

[0067] The glycerol solution which the marker particle mixed in drawing 7 is 
installed on a rotating disc, and the result of having computed the movement 



magnitude of each particle by having carried out the image processing of these 
60 particle track images for the particle track image which photoed false flow to 
drawing 8 by the particle track pursuing method is shown. 
[0068] 

[Example] Assessment of the engine performance of the floating instrumentation 
system of the fluid concerning this invention which performs the image 
processing by the particle track pursuing method to below, and the applicability 
to the real floating field of a floating instrumentation system are explained. 
[0069] 

[Example 1] (Assay of movement magnitude with a jogging base) the true value 
to which move the particle image drawn on the plate with a jogging base 
(mechanical stage), and measured the movement magnitude using the 
70mmx70mm calibration plate, and it was made to move actually, and the 
floating instrumentation system of this invention — Streak - the measured value 
measured by law was compared, the error was surveyed, and the precision was 
authorized. 

[0070] What has the lattice point of 64 points was used for the calibration plate 
into the 70mmx70mm rectangle. Therefore, the measurement range is 



70mmx70mm. The image pick-up scale factor at this time was 5.01 1 mm/pixel. 
[0071] The particle used three kinds of particle diameter 3 and 6 and 8pixels. A 
particle image picturizes the particle by which seeding was carried out into the 
fluid with laser sheet lighting, prints this out by the printer, and is made to paste it 
up on a plate. 

[0072] The image pick-up of a particle image picturized the particle image of 
plate Kami installed in the criteria location (criteria location in a calibration), and 
the particle image after migration with a jogging base with the CCD camera, and 
made this one image pair. It is equivalent to the approach which depends 
calculation of a movement magnitude vector on a cross-correlation method from 
an image pair, and is used for measurement of a real floating field. The number 
of vectors computed from the image pair of a lot was about 150 pieces. 
[0073] Migration of the plate with which the particle image was drawn was 
performed in the x directions to every 50 micrometers 1 0mm, the image of 200 
sheets was acquired, and 100 sets of image pairs were obtained. A 
cross-correlation method and normal-distribution approximation (Gaussian 
peak-fit) performed movement magnitude calculation. 

[0074] Since the diameter of an optical fiber single fiber of the image guide 26 



used as the minimum spatial resolving power is df=3.12pixels, it becomes d 
tau/df =0.96 and 1.92 or 2.56 times from the frequency analysis result of an 
image guide 26 to three kinds of particle **** 3 and 6, and 8pixels, respectively. 
[0075] drawing 9 , drawing 10 , and drawing 11 — respectively — Streak — the 
bias error (Bias error), the dispersion error (RMS error), and probability density 
distribution (PDF) of the measurement result by law (the particle track pursuing 
method) are shown. In each drawing, b and the curve of d tau/df =2.56 are set 
[ the curve of d tau/df =0.96 ] to c for a and the curve of d tau/df =1 .92, and the 
axis of abscissa in drawing shows true movement magnitude (Actual 
displacement) [pix.]. True movement magnitude can set PDF shown in drawing 
11 to 0.5pixels(es), 2.5pixels, and 4.5pixels, respectively. 

[0076] Drawing 9 shows the bias error of the measurement value of the amount 
of x directional movements when moving a plate in the x directions, and can 
move aside from the number of migration pixels of a true value, it is what plotted 
whether **** shifted and measured, and the bias error shows the difference of 
the average of all the vectors of measurement within the limits, and a true value. 
[0077] moreover, drawing 10 the root mean square of the measurement value 
of measurement within the limits - taking — drawing 11 — true movement 



magnitude - the time of 0.5 pixels from a left crest, 2.5 pixels, and 4.5 pixels — 
Streak - when the measured value by law (the particle pursuing method - TV) 
set the whole number of errors to 10 (%), it be show whether actual 
measurement be carried out in the place which shifted from how many pieces 
and a true value (%). It expresses that there are few errors as the width of face 
of this crest is narrow. For example, having measured with 2 pixels means about 
nine of 100 pieces at the time of 2.5-pixel migration of the crest of middle. 
[0078] As shown in drawing 9 thru/or drawing 1 1 , by the particle track pursuing 
method, the difference of the measurement value by particle diameter seldom 
appeared, but bias error beta was settled in the range of abbreviation 
beta=**0.2pixels. 

[0079] Since there is a clad (permeability is comparatively low) which does not 
let light pass in the optical fiber 25 used in the experiment, lack of a particle 
image becomes a cause and the oscillating effectiveness is seen. It is because 
the value change which depends on particle **** is not looked at by the bias error 
shown in drawing 9 but sufficient information can give this to it to the diameter of 
an optical fiber by the particle track. 

[0080] Streak shown in drawing 10 ~ although the dispersion error by law did not 



depend on particle ****, either but the almost same value was shown, it became 
a comparatively big value. The configuration of a histogram not serving as 
normal distribution, but becoming horizontally near distribution near a true value 
results so that it may understand, even if this looks at probability density 
distribution. This is because correlation of a particle track image is taken, namely, 
the brightness value of a particle track is comparatively gently-sloping, and is 
because correlation coefficient distribution of the image will also become 
gently-sloping. By the particle track pursuing method, having taken the large 
particle compared with the usual particle image enlarges dispersion error 
epsilonk (0.3 - O.Spixels). However, not carrying out the method of generating 
with error that PDF roughly breaks by particle **** was shown, and it was 
checked that the precision of a bias error is effective. 
[0081] 

[Example 2] (Measurement by the real floating field) The floating field 14 of water 
(fluid) was established in the experiment container 15 as shown in drawing 1 , 
the axial symmetry jet was measured with the floating instrumentation system 10 
of this invention, and the engine performance of the floating instrumentation 
system of the fluid in real floating field measurement was examined. 



[0082] Examination of the floating instrumentation system of a fluid was 
performed as compared with measurement by the Doppler type laser 
velocimeter (LDV). 

[0083] Containers 15 are 250mm long, 250mm wide, and a cube container with 
a height of 250mm, and gave the rectangle jet outlet (nozzle) 16 in the center. 
The 3-30-micrometer globular form polyethylene particle was mixed in the 
working fluid as a marker particle at this using tap water. As shown in drawing 1 , 
once accumulating a working fluid in the lower storage tank 50, by the 
experiment facility, it was supplied to the upside reservoir tank 52 with the pump 
51. When the working fluid supplied to the reservoir tank 52 makes a part 
overflow through amount of batches 52a, water head pressure (head) is kept 
constant, a working fluid should pass the honeycomb and the mesh (not shown) 
the floating adjustable valve 53 and for rectification from the reservoir tank 52 - 
rectification-ized, it was made to spout in the facility container 15 from the axial 
symmetry nozzle 16, and the floating field of water (fluid) was constituted. The 
water overflowed within the reservoir tank 52 is returned in the storage tank 50 
through a return pipe 54. 

[0084] On the other hand, the contracted vein of the surface ratio was carried out 



to about 10:1, the stable velocity distribution was acquired, and the nozzle 16 in 
a container 15 conducted the floating measurement experiment of a fluid(water) 
in the condition of becoming a stationary mostly. The system of coordinates of 
the floating field of water (fluid) made the core of nozzle 16 outlet the zero, and Z 
and radial were set to r and they set [ the direction of a jet (the height direction) ] 
the velocity compornent of theta and each direction to w, u, and v for the hoop 
direction. 

[0085] Carrying out outlet velocity of the working fluid from a nozzle 16 in w= 
165mm/second, the number of lei nozzles to the outlet diameter (D= 7mm) at 
this time was Re=1 155. The water which overflows a container 15 was also once 
returned to the storage tank 50 through the return water piping 57 after the 
carrier beam by the receiver tank 56, and constituted the circulatory system 58 of 
water (fluid). 

[0086] The probe 40 for visualization as shown in drawing 4 which the scan 
optical system 13 as shown in drawing 4 , and the picture transmission means 
22 unified was used for the scan optical system 13 and the picture transmission 
means 22. As shown in drawing 4 , the probe 40 for visualization was inserted in 
the floating field 14 of water (fluid) from the upper part of a container 15, and the 



laser beam was supplied to the vertical from the upper part, and photoed the 
particle track image by the fiberscope 44 from the probe 40 which counters. The 
distance of the laser sheet 17 and the light-receiving side of the probe 40 which 
counters is about 30mm. The measurement range is rectangle range which is 
range where the lattice point of a calibration plate exists and which is 
10mmx10mm, and the thickness of the laser sheet 17 is 2mm. 
[0087] To the diameter D of a nozzle 16, the potential core of the jet of water 
disappeared, the photography object domain was called developed region, and 
the momentum exchange with a perimeter fluid set it as the field of Z=6D to 
which a comparatively gently-sloping velocity distribution is acquired even in the 
jet shaft center. 

[0088] The picturized particle track image performed the image processing by 
the particle track pursuing method (Streak law) with the image-processing 
means 35. 

[0089] On the other hand, measurement by LDV (Doppler type laser 
velocimeter) also went separately in the cross section vertical to the jet shaft in 
Z=6D. These measurement assumed that a floating field was regular, adjusted 
the flow rate to the same conditions, and performed it. 



[0090] The jet outlet conditions (Reynolds number: Re=1 155) measured by 
drawing 12 by LDV are shown. The curves d and e of mean velocity (LDV 
(mean)) and fluctuation velocity (LDV (rms)) are shown in drawing 12 . Moreover, 
the velocity-vector map near height direction Z=6D is shown in drawing 13 from 
nozzle opening by the particle track pursuing method. 
[0091] By the particle track image, the rate was computed by the 
above-mentioned particle track pursuing method (Streak law). The 
shaft-orientations mean velocity distribution f in Z=6D and its fluctuation velocity 
distribution g are shown in drawing 14 and drawing 15 , respectively. An axis of 
abscissa is radial and axes of ordinate are shaft orientations. The measurement 
and LDV measurement by this system performed the comparison of the result 
which extracts this field, performs statistics processing and is depended on two 
mensuration from the shaft-orientations mean velocity distribution h which is an 
acquired velocity distribution, and its fluctuation velocity distribution i. Being 
computed corresponding to each particle, and since it is seeding of low number 
density comparatively, the locus image by the fiberscope in the floating 
instrumentation system of this fluid has little rate data obtained from the image 
pair of a lot. In order to compensate the little of such a statistic, a lot of images 



needed to be photographed in the direction of time amount, and the statistic was 
computed from the image of 16000 sheets in this experiment. The 
shaft-orientations average and the fluctuation velocity distributions f and g which 
are shown in drawing 14 and drawing 15 were statistics 1000 [ about ] on each 
point. In measurement by LDV, an average and fluctuation velocity data were 
computed from the statistic 3000. 

[0092] if the shaft-orientations fluctuation velocity distributions g and i shown in 
drawing 15 are compared - Streak - LDV and the inclination of law 
corresponded and bias was looked at by the magnitude. About this amount of 
bias, the larger amount than what was shown in the experiment by the feed bar 
Cogging base) was computed to LDV at the jet core (r= 0). Having the 
measurement value which can be trusted was checked about the 
shaft-orientations mean velocity distribution f and k shown in drawing 14 . 
[0093] if a fluctuation velocity distribution is referred to - Streak - law - LDV and 
an inclination - the fluctuation velocity distribution i of LDV of a match - 
comparing - Streak - the fluctuation velocity distribution g of law showed the big 
value, from the irradiation time of laser - converting - Streak - law was set to 
0.63pixels(es) to LDV in r=-2.2mm, it has fitted in the precision by the feed bar 



assay result, and it was confirmed that it is comparatively hard to be influenced 
of shear. Near a jet shaft, fluctuation is large and it is thought that minute 
fluctuation in the small field of a velocity gradient was not caught. Moreover, 
although fluctuation shows the big value compared with the perimeter in r= 
3.3mm and near r=-3.3mm, mean velocity is a low speed and it is thought that 
fluctuation is measured greatly relatively. 

[0094] drawing 16 and drawing 17 — Streak — the shaft-orientations rate 
probability density distribution (PDF) j and k of the measured value in the jet 
shaft center (r= 0mm) and shear layer (r= 2mm) by law and LDV, and I and m are 
shown, respectively. Both mensuration is known by that measurement is 
performed normally in a jet shaft center. Moreover, it turns out that measurement 
is mostly performed normally also in a shear layer (r= 2mm). 
[0095] The applicability to the engine performance and real floating field of a 
floating instrumentation system of this invention which perform the image 
processing by the particle track pursuing method is checked, and floating 
distribution (the rate of flow and flow direction) of a fluid can be measured with a 
sufficient precision to accuracy with this floating instrumentation system. 
Thereby, the floating condition of the fluid (coolant) of the floating field in the 



severe closed space which differs in the external world and the environment like 
[ in a reactor pressure vessel ] can be measured to accuracy. 
[0096] For example, the floating condition (the rate of flow of a coolant, a flow 
direction, a reactor core flow rate, stagnation condition) of the fluid in the 
downcomer section of a reactor pressure vessel or a reactor core shroud 
(coolant) is measurable. In a reactor pressure vessel and a steam generator, 
since the heat floating conditions in the heat exchanger of a thermal power 
station plant etc., i.e., flow, can measure to accuracy heat floating (floating 
distribution) which is not uniform with a sufficient precision, the thermal fatigue of 
a component part and the progress situation of corrosion are manageable to 
accuracy. 

[0097] Since the thermal fatigue of a component part and the progress situation 
of corrosion are manageable by specifically measuring the stagnation condition 
of the support plate supporting the heat exhanger tube of the interior, such as a 
steam generator and a heat exchanger, and the fluid formed in a gap with walls, 
such as a heat exchanger, etc., it can operate efficiently safely more. Moreover, 
since the floating condition of a local fluid is clarified, the design of optimization, 
such as a heat exchanger, is attained. For example, according to the floating 



instrumentation system shown in drawing 5 , also in a complicated gap, it is 
easily measurable. 

[0098] Moreover, for example, in the floating instrumentation system shown in 
drawing 5 , the laser sheet of the 2-way which irradiates and goes direct so that 
a sheet-like laser beam may be intersected perpendicularly with a 2-way from 
the scan optical system 13 for laser sheet formation is formed. By picturizing the 
pathline of the particle on two laser sheets which intersect perpendicularly with 
CCD camera 24 through each picture transmission means, a motion of a particle 
can be measured in three dimension, as a result the floating condition of a fluid 
can be measured in three dimension. For example, what is necessary is to install 
the body casing 45 (to refer to drawing 5 ) in the downcomer section free [ the 
migration to either / at least / radial or a hoop direction ], to measure the 
movement magnitude of each particle to accuracy in three dimensions, moving 
the body casing 45 in radial or the hoop direction of X, i.e., the direction, and the 
direction of Y, and just to measure the floating condition of a fluid, in order to 
measure the coolant flow rate of the downcomer section of a reactor pressure 
vessel to accuracy. 

[0099] In the floating instrumentation system 10 shown in drawing 1 , by unifying 



scan optical system and a picture transmission system, as shown in drawing 4 
and drawing 5 , it can apply as a flight model and can apply still more broadly. 
For example, while being able to attain high-degree-of-accuracy-ization of a 
floating calibration of an internal pump by measuring the stagnation condition of 
the circumference of the internal pump in a reactor etc., the thermal fatigue of a 
component part and the progress situation of corrosion by the stagnation of a 
fluid etc. are manageable to accuracy. Moreover, it is also possible to be able to 
apply to the flow test measurement accompanying enlargement of the structure 
in a furnace accompanying amelioration of ABWR, the formation of shroud loess, 
etc., and to contribute to rationalization of a design and the improvement in 
effectiveness of a plant by high-degree-of-accuracy-ization of the further 
information. 

[0100] Moreover, this floating instrumentation system is applicable to 
measurement of the floating field of the fluid [ structured division / in contact with 
a stream or an air current ] of the circumference of equipment and materials in a 
marine vessel, the aircraft, or an automobile. In this case, the floating field may 
be formed with the gestalt opened in the external world. 
[0101] For example, the floating condition (the rate of flow, a flow direction, a 



flow rate, stagnation condition) of the stream in the pars basilaris ossis occipitalis 
and wall surface of a marine vessel is measurable. Measurement of the floating 
condition which is not uniform can also be performed with a sufficient precision 
to accuracy. Progress situations, such as adhesion contamination of the shellfish 
which this generates in the stagnation part of the flow in a pars basilaris ossis 
occipitalis, a wall surface, etc. of a marine vessel etc., and fatigue of equipment 
and materials, are manageable to accuracy. Moreover, it is also possible to 
attain optimization of a marine vessel design and rationalization by flow test 
measurement accompanying the configuration of a hull, structure, a member, etc. 
[0102] Also in the aircraft or an automobile, it can contribute to the optimization 
design of improvement in the speed, transit safety, etc. by flow test 
measurement. 
[0103] 

[Effect of the Invention] According to the floating instrumentation system and its 
measurement approach of the fluid concerning this invention Visualize the 
motion of each particle distributed in a fluid with a laser sheet, and the particle 
track image on the visualized laser sheet is picturized with an image image 
pick-up means. Each rate of flow and flow direction of a particle in a complicated 



floating field can be measured, and floating conditions, such as the rate of flow of 
a fluid and a flow direction, can be measured with a sufficient precision to 
accuracy because compare and analyze the particle track image of 2 picturized 
time of day (minute time interval) and it carries out an image processing with an 
image-processing means. 

[0104] Since measurement of the floating condition of the fluid in a complicated 
floating field can be performed by processing a particle track image with an 
image-processing means according to this invention, the floating condition of a 
fluid can be measured from a remote place, and the floating condition of the fluid 
in a closed space severe in environment can be measured with a sufficient 
precision to accuracy. 



[Translation done.] 
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[0 0 16] *%nic{|«SIE^OKItrtN-ai^tt. ±2 
10 illci/-h«icMLTU-fv'-b^ML, CO 

9m lxk7t4' m$mm<D&m> w-y 
ftttZMtb, c<Dyt77-<'y®&m&<Dm&'Vi-y* 

mm®icmzwi&<Dffim*ftm?5J5mx'$>z>o 

[0 0 17] 

7-L.&£zs?(Dtmm&v-nffiftm^^xmttmn 
*pmLxmmt%o 

lo o 1 8] m i«, *«WK«*at(*oaiaBrtH-i»5'^ 

[0 0 1 9] .*»WO«fflc©aaBifrflJ'>^f-A 1 Ot 

-essm* u— •fftfcsig-r w if*/ 

30 £]i$ft8 0 8 nraO^ftU- W3tt«t LTfflV^ 
[0 0 2 0] U-+f'fgSgH 1 1 ft»5«a*tlfcb-1f 

fiiSIl l*6«fi*nfcW-1f3^illEtt:0«IWIl 

[0021] M^o^Kit 1 4tt, nntwi&z&Kt 
40 j:^..aawii 4«> ?s<*co9-5^ / >*< tt,^ 

»W*OJ:5ft, 7]<«[-»f>Stfflji:ft«f SWigSPffi-^att 
g t» 0 tc fett § t 4: t § c i: *<t? t 

§o **3, 0 l fcw-rJKaHTW:, ffiftOdKHM 1 4 4 

50 7j<(D^iljil 1 4*MLTV^. 



5 

[0 0 2 2] jtSft^l 3(4, ffi.ft<Dti5M®l 4lCl> 

uTaf*aawii 4h:u— tf h 1 7*«*tTw 

£>„ COU- ■fS'-h l 7fcJ:0SS#OilSIWil 4^pJ 
flMfcl/Ci^*. 8itiJfta!l5/*f A l 0T*(4, ffitttiffiftil 

i 4©*ff**iaioami»»jii*«ffr*fci6tc, »& 

[0 0 2 3] jfegft^&l 3(4, B 2 fcaVfJ: 3 U- 

7^1 2fr5**ai£tf9ft*#oTtli*t3ft*. til 

[0 0 2 4] COftft, 3(C(4, gftffl7 

7^ 1 2fr6lti»'*nfc^- J !ftr-A^5>b— tf->- 

LTvuyxmi s^ffl^ens. 8 
(4, ^fe-A^MtTi^-ife-Atcnu^-b 

£-e3 2fe©3'J^-hUyXl 9 a, I9bk, 3'J 

— tf e- u— tf j/- f 1 7 ©j?**iq]-e*?j i ~ 2 m 

m@|g{C7XXh^^g^b>X2 0 £, U— 9*5^ 
-h i 7^fig-TS>'U>'KU*;l/UyX2 1 fcfcffl* 

U^-hUVXl9a, 1 9 bfcti. flfctf, 

[0 0 2 5] j£4tt*Sl 3*^->-h^tcra#t^n^) 

u— vit\c «t o , «tf*©&£»j« 1 4 tc u— ys/- f 1 7 
jwe«?tu .stftowirti 1 4 2^>- htttcojtifk^n 

IT, U— »f->-F 1 7©fflJ£«Hl 7 afcftfo 

-f 1 7 ±omimmmmmm%2 9-?mmmt%:z> 

£9(c, j8fflflnS;tf7*2 3tfSaflSftTVSo 

[0026] mmmm^® 2 9 u— >f ~>- f 1 7 ± 
©Bfg^iSM-r § i5«eas#a 22^,2 ^s^hisif 

*rr*. H«eai#a2 2«u-ifs/-n 7 ±02* 

ea«-&*J;9h:ftoTV'»«. BiMs3M#IS2 2(4, 0 

3(c^-r<fe5(c, aw5sro*©)t77^2 5% 

6ffJB^&ftS, -f^-v^V F2 6(4, 0 3£jjVr«J: 
5(C ^fe5»n/c7 l c7 7-r^2 5 ©PWBfl5©&{ftffltfIE 

5S(c wis-r s «fc 5 (cms ? nr *fe ?) n« o -r * - 

^ F2 6©B»4¥ffilctt±(fe>ft, WhXsVXZ 7 

(c«to-^<D7 7^/^®(cM : ^nrciii^#iffi^7 

CDA>t72 4©J©7 7^/^E$T?fcj2ILT^<5 6 7 



(4) . ^ 2003- 8 4005 

6 

r^/^ffi^e>*J?«v 1x^X2 8^STCCDA^5(c 

1"*(C, L/-^>- F 1 7 ±© 2 #jcS7#l8frBll£iig 

[0 0 2 7] -^-i^f-f K2 4©B§5©Plf? GK£ 

oTft.JO, fi©ji;i7j(4, F2 6(c^fe 

e>nsft77^2 SOl&MSfcM^fcJ^TBli: 
£o ft77^M2 5<037g*0'M*3t»O3t»fl5 
10 (4, ft77^2 S^fitoSHK^ftSn, #'37W 
®rtT-(4BJ33£(4Hi£:&5o LfctfoT, 
-YF2 6(4ft 77^2 5 ©378*** (tf^-t;l/) 
fctT.©3SniKM!ll60«'J>*fflfcft*. 
[0 0 2 8] ^fiSHM F 2 6 (4, 

OflHEfcEfl* ntV*. F 2 6 T*ft 

So 

20 [0 0 2 9] Sfc, ^77^2 5RBfc*7-y Fffifttf 
Sfcf £nT^5^->W F2 6T'(4, ?v-y F©ff 
^^SK&Sfcx 100^77^^2 5<037(CAofc 

tiuoWSfbWfr* ^©3yKvXF*Mg"Ft-So c«fc 
fe, ^)«-^VK2 6(4, ftafflBM^J«3fi*8fli L 
fefc*/ FOJP**Wiar=8:^ci:^*ffi:* 

[0 0 3 0] -7^, SW*<OSHftlf-MS/^'ri» 1 0(c(4, 
^S^ayhD- J^St3 O^fiixSnTV^o ^ 
30 -T^ y^3yFu-;V^3 0(4^-<5>'^'X^v f aL- 
5 3 1 fc^V^a-r-f+Fs 2fc^e>^fig^n§ 0 

i/y^Dt^f 3 2T[p]$£-<h-=>TC CD*^72 4^ 

*s 3 o (4, v-4?mmm. 1 1 tswwsfc lx<d c 

CD*^72 2il(D[B]^toT|g®J$^?.^d(C^o 

[0 0 3 1] CCD*^5(C(4, Px(4, ifi*^*RfSW 
40 lcmM*ft-Dffl&&V GA (640X480pixel 
s) , 71x-Ab-F30Hz, 8t:-yh ftyyn2 
5 6^11) ©feO^ffl^^nSo CCD7J^5 2 4(4, 

u-<fis-h 1 7 ±o 2 &fcm?%mw®% \s-*fmsi 

SHI 1^6©»«*-f5>^fcHJW*fc-DT»«LT 
«30, C CDA^72 4T'ffimZtirc7i-a7<DWm{E 

5tt?<W»ia»ffifc * 0Mi*©Mii^ l 4 (cfe(7i.^ L {* 

©«Eai^^n73io)ATOJ^nSe ccd*^7"2 4 

(c(4, I?iJk(4, a^tta^FS (7a^U>y^77+^ 

50 y) <D**?&m\i^tu mnsntzc cd©^^-t 



(5) 
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[oo3 2] ztc. mm&m^®3 5ti\ CCD£^ 

2 4 % <D 7 -f U VWrnim*? is 5 Mtt & A / D g 
MstLTfDyV-LVU^-YZ 6 C(D7ls- 
L??/^- F 3 6fr6x-7#;MItflfI^#P C I MX 

^©m^3 7^LTiMe»n, x^nsMinirs 

3>tfa-^3 8 fc^Wf^o 

[0 0 3 3] C(Dffi.&<Dffi.mVMi' XrA.l Ofi, C G 10 
D*^72 4tT7tD^U HMLS^KS 5fc 

[0 0 3 4] 01 (c^-T^O&itijf+iJv'X-fA l 0(i 
tfB 5 ic*? J: 5 <f >>- bffl^Sft^ 1 3 £ 

w®mm&2 9 »{k^n^.o 0 4 . 
& 1 3 t-MtLrcwmmm®2 9<Dmm&m^mz 
2 mm i mm.mmm^u--f 4 0 1 ur^e^ti 
T^3 0 7n-7 4 oft^^B^n/c^^ (Ji^ffl 20 

SWf^*) 4 1^5 4 5° B^T^Xi^ 3£«fcD7tB 
^ftlfi®fie^#g%m^t/c7rf'M-Xn-7 p 4 4 
(CTBfMSj^nSo El 4 (A) CO-fWtUfcilK 
»08MWtMS/Xr j**H l ?>S2§ l 5 

t fc^s^^^ 1 3 1 mmmm^m 2 9 <om»&&£ 

@2 2 k**flcir— >y^4 5te««LT^«. 

[0035] wfroaaBirai^Fffisittw-r^.- 

[0 0 3 6] 01 tC^^nfc^cD^KH-Syt/XxA 1 
CCDJ3*?2 4lCl&?<Dtfim> (flBflO 

s 2 ^jos?«i»ii«%iWfi^^-^i:**LTia« 

iiSfrffi (Particle streak track 
ing velocimetry) t: n?^ 0 C Oifi?®! 

[0037] n^mmmt. u-if->- h 1 7±<d?h 
mm wcwmu tm#iB*&to*0M¥ (wmm 

&$Mkztirm&*&'DV*'bJi' (mm) ia^ti 

&a^w*Hfc035H£*nfcj©Bi!fitti 7 a 



ftffl 2 0 0 3 - 8 4 0 0 5 
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[0 0 3 8] tt?«lttBftfr5zk (»tt*) 

[0 0 3 9] LfrU 4*—Vi3'< K2 6£5ILTC C 

^©ffS/<*-:/7£W-Tr!S:<, >T K2 6<Dft 
7WM2 5<DBE?"J («*#©$§£) <0»fi/^->t. 

[0 0 4 0] *-<lT\ *«Wfcff*«M*:0«SU»IHIB3/^ 

[0 0 4 1] 1SJS2 0 0 0-2 0 7 4 7 6*f&fgfc&3 
HTOjiltSiJ (P I V) }S*3SfflLfct><DT»a60^ C(D 

p i v&T-&, aea^ffl^enso jimop rvs 

T*^ WfVmtM ^-^K2 6C«7r^2 5 
[0 0 4 2] Sfc> «/J>oa»lffliMllli6*4A*3t7 r>T 

^ 2 5 ©*s»t-*«fc o fes^fts^^t < -r § i:s e 

gffil? t iBfffcS: 6 - 8 11 %> % C i: * L 1/ ^ ^ n 5 
46, #BBSSli»*JtBWl*#<9l6*<Ttt46a<a 

[0043] c<Drc&b, mwm®*m^z'p i vst- 

tc «fc t) , iSfiLhKS^fti*** <ii^ci:^<, ^ 
•3g/J^r B % v ilfi^^S7t7T'i'/^2 5©#W;M 

Lr+^affi^ofliw^^sc ttfTirs. cos 

^il^iB^j* (Particle streak tr 
acking velocimetry) 



(6) 
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[0 0 4 4 ] *jffivm^*&?wifiesffi"ett.* ^<r> 
co o 4 5] ( i ) j^xomk mmmm) 

^^-y*VK2 6£iILTCCD£^2 4Tlt$ 

?o 10 

[0046] 9t7 7j'mw<D&m'W-yffi£<Dwm 

[ o o 4 7 ] mizt- ? <D%iL^-mm\m xttzt 

*, e^niafflfiiixfi, 

[fti] 

x= 8 + n ( 1 ) 

■es^rnso 20 

[0 0 4 8] T ftfc>%, jfliJ^x-^O'lS^flSfWixtC 
f4, (tt7r^E7!l<0»fi><*-l') CK 

/w-ytt&umtt. sNtt*^<, imm^vm 

[0 0 4 9] Sl*l^3i»>JBJlS*ft1»LTa*r*i:x 



[0 0 5 0] 30 
■[»2] 

x' (k) = s 1 (k) + n« (k) (2) 

i=i*i®e?i*. 



k : itfegOta^aiftWcDkSgCVf^-fe^tt*^. 

[0 0 5 1] «a?«l»H«Oll»fciBV'«T, BlDfcA/JS 

fludfrrainMWSx (k) 3 

[0 0 5 2] 40 
[ft3] 

[0 0 5 3] 3t7r^r^E5lJO»fi-'^-^s{±, II)- 
©jbO«C*DM{g(Z>ffgfiS£ft!K iin§¥Wl> (M 50 



<&m 2 0 0 3 - 8 4 0 0 5 
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7r^ <E5tJ0»*> nfli i: ft * . 

[0054] cfD^msaatcfcoT, treats 
«Wfit4iffiffii^ft:tt>feors{^-tftsor\ fflBB 
tm[«iafc» tTWSft/e^-eft < , zm&maa^T 

^W^ft¥Skfti.o 

[0 0 5 5] ( 2 ) -imt • B$f biiftfta 

T, -ffl<fcfc*tf*SWtH«afiS*fT$. ?K (SB*)' * 

OR****, »KWftt»©fc«ft&1\ Wfil«ftli« 
ilLTlBit^n^ (0 6 (a) #SSJ 0 cotw&wta? 
ffi^rfi|{kMStTt, ^77 Kte«fcoTtiJIB©5fcfc 

nfe*$H*iHi*snft^ (06 (b) » . 

^TOlijR* 0 35)^2 5 5fCSSlt:St>OT*S (83 

* 2 niaig-r c k T*»»f $ tiitmstm^ mmt, mm 

SB) o 

[0056] ( 3 ) >7"fc<fc :rfe?H»ffl«o* 

SftofcSSS^fcligijOS^^oltSo 

;L>?£ (Particle image centroi 
[0 0 5 7] Xo Yo U 

[»4] 

( 4 ) 

[0 0 5 8] (4) ffiHttfffc££BSM$ffift0ffti] 
(1) ~. (3) ott?H«ia«yifcJ:-3T> li?»ia 

[0 0 5 9] ft*3, Sft{*<D«t»)ffS^*xi» 1 0©C C 



(7) <&M 2003-8400 5 

11 12 



#«• 5 2 ^Tc^fiM^^r^^^a l fee? tug 

[oo6o] mnmrnm t * sia? o»»*«i««, * 
/?(^',y',i,»7) 



[0 0 6 1] 
[»5] 



[0 0 6 2] 



n-l m-i 



* « [»6] 

1 n-l 



I > ' 2 

n-m nm 



V) 



(5) 



•(6) 



[0 0 6 3] S5**tf£6fcjjVrJ:3fc, WiR 
4:*. fflH«R#l0fc*fc2Htt4£±K:HRU ffi 
[0 0 6 4] ^5fc«fctfS6K:*5^T, I i i 

AX' , AY' IcmmtZa «BtO^?nXm 

[0 0 6 5] (5) VfVtizMmmQmti 

(4) offlBBfWtJit), tfc?iS«®$ll!i*<Dfi»ffitf 

^^frifiiW (Gaussian peak— fit) tc<t 
[0 0 6 6] (6) $51 

W±<D? life «t K) » t^LttBftQ&Hfttf t? * 
(pixel) «IEfc:<fct>^£ISlT' 

[0067] El7tch ^-tt^ioax brc ^ U -tr U 
y»»*lHllSP3ffi±tcfftlBLx JJWttW*Mn«JMgLfe 
*fr?IM*H«*, 0 8 C OtS? «lWiS« 6 07WA 

[0 0 6 8] 



20 fc c ttf»ltiiJv / XrA«^»l-\©jgffl14fco^T 
[0 0 6 9] 

T^ftS-S* ^(D&WlWSn Ommx 7 0mm£>$$IE1£ 
MOBllHB^f-ixtei'J Strea kffitc J: OfHULfe 

30 [ 0 0 7' 0 ] fSiHS^x 7 0 mm x 7 0 mmOj©Krtt 
6 4j£©«?Mi*£*St>©*ffil/\fc. LfcAbt, It 
$JiBH&7 0mmX7 0mmT?S5, 'CQfctQ&Mfeffi 
${±, 5.011 mm/p i x e 1 T'&ofco 
[0 0 7 1 ] te^titft^Sa s 6, 8pixels£D3 

[0072] tt.=pmommt. mm 

[0 0 7 3] S?«©ffl3b^nfc¥SO»«I«. x^C 
5 0 jimOO l 0mm$T*!f</\ 2 0 0 &<DlHM:lx1# 

i ooffioiMW**i#fc. ^®iawtB«, ffistira. 

S^oi^IE^^^jfiftl (Gaussian peak- 
50 fit) tCj;0tfofc o 
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[0 0 7 4] gtfr<Dm&#Mte£1&&<( F 2 

a*?«r*S*^6df =3. 12pixelst^O 
T\ 3fi£l<Df5^<iS3, 6, 8 p i x e 1 stJtL 
T, ^n^tldr/df =0. 9 6, 1. 9 2, 2. 5 
6{gil&;f>o 

[0 0 7 5] 09, 0 1 OfcJctfB 1 1 K^ft^fft, S 
tre akl CK? «iSBi»ffi) te**tH»B«Ofi.9 
if (B i a s error), f£P>OtHM (RMS 

error) feJctfttspfSa^ft (PDF).*St. 
&0{c;f3V">T, di/df = 0. 9 6©ft*fc£a % d t 
/df = l. 9 2©«^b, dr/df = 2. 56<0 
ft»*ci:U 0^^ii&K<D&l» (Actual 

displacement) [p ix. ] *7r;"f 0 0 

i nc^$nfcPDF«, n<D®mmti i zft : eno . 5 

pixels, 2. 5 p i x e 1 s , 4. 5pixel 
stents tOTSSo 

[0 0 7 6] 09te, ¥ffi*x£lfcl«:gfl&£-frfci:*© 
x7aln]^»J»«ltWIcD{aO^M^*Lr*5D, Kfl© 

[0 0 7 7] £fc, 01 OttfHUBHrtOttilfilO;!'— 
h fcofefeO^feDv 01 Hi, n<0®W}M 
t>\ &©iJj*^0. 5e^t;K 2. 5tf^-fe;V, 4. 
5 e^-tr;KDB#, StreakS dKPjBMftt : P T 
V) fcJ:S«SM#, ±f*©8SM»*l 0 (%) tLtc 
B#, fafS (%) , KfflJ:0mfei:c5-e*ISWai*ti 

5 fcf*-fc/M£«j©B#, 2tf^-b;bi:ft!S'JLfc©(i, 10 
0{B© 3 "51519 IB fci^S C fcSriftiftf So 
[0 0 7 8] 09&^L01 HC^£tl/c<fc5tC, gff 

hi-, m*)mgpimp=±o. 2 P i x e 1 s<D$ga 

[0 0 7 9] ^T'ffl^^7 7^^2 SfctiftfcilS 
^•o^ff^lKHfcftO, JHW»*tfJI5nS. 09(c 

5. c t tfX* t 3 C t J: 3 0 
[0 0 8 0] 0 1 Olc^tiTzS t r e a k 

5 o * mm ftLTMmcfcz -rmm-orn^ u t ^ 

T£9, -ffcta^ tt^*l»<D»ftffltttt*fl^&J£5>fr 
T $ t> , * ©HfcOtB Hftft^ & £ ft & f£ 5> & t <D 



3) #§32 00 3-8 4 0 0 5 
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iie k £;*:t<f.5 (0. 3~0, 5pixel 

s ) c t , Kmft&ic <t o t pdf < fijn-c 

[0 0 8 1] 

[usstfij 2 ] (saawi-eottao 0 itit * 3 

^S§§ 1 5 P^tCTK (AH*) Oil* 1 4 *RttT*J«fl 
[0 0 8 2] 8S#0«tE«l9ra!l'>X-fAOttf*tt^ >*y7 

73v- tfaianH- (ldv) fc^stHifctHSL-cfi^ 

[0 0 8 3] nmi 5tt, 0Rt£SK2 5 0mm, ^2 5 
0mm, iS£ 2 5 0mm©£7J{*^3gT'&D, 

^ps^tBn (yxw 1 6 £if fc-e/co fHttatftta* 

Jl7K^fflt/\ iintChU— 9— &?£:LT3~3 0 (im 
20 O^^Ux^l^yli^^MAL/co 0 nc/fs-fefc? 
(c, MMBHM-p(d:fmSKfM&. T0|5Olf?@^y^ 5 01c 

5 2'\{&S£Lfc„ U+f'-/^^^5 2JC^$nrcft») 
ilM^tiU^-^^V^ 5 2A^MS]slli5/Vl/7*5 3, 

30 2frx*-^-7u-Lrz7mmK>W5 A?zm.x%rmz 

5 o\HicmZti%° 

[0 0 8 4] -73, 1 5rt<D/X;l/ 1 6tiffittJ±^ 
(?SS(C^5mi^*3V>T^ (7j<) 0«JilthiB!j*»% 

tTofcc 7K cat*) <ommm<ommMtyx^\ 6tup 

1*1* r, JSlTf^^e, ^-tl^tlO^lPJOiSa^Sw, 
u , v t L/Co 

[0085] /X;n e^eof^ij^ftottjpJis^w 

40 = 1 6 5mm/&^U ElOfctOfflPiBS.CD- 7m 
m) £#-f£W/X;H&&R e= 1 1 5 5-efeoTCc 
gg§l 5**- /V-7D-1-§7j<£>, ^t^y7 5 6T 
— SSttfcft, RD7KSSW5 7*gTlrtig^^^ 5 Otc 
K£ft, 7j< (»[#}' <D«SS5g5 8*«S«Lfeo 
[0 0 8 6] jt*Jt¥3R 1 3 fciS^feji^S 2 2 tt, 0 

4 tc^-r j: 3 4j£a^» 1 3 2 2 1 ^ 

-i*{kL/c04tc^-r«fc-5*Rr^{t:fflXn-7 4 o*ffl 

t^fco 04 t^-r^^tc, nrmfflXa-X4 o*si§ 

50 tf3fctt±*fr6»iIfcfe\StU MW57D-74 0 
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?j£##£-f ^fEBT'S^ iAIflOmmxiOmm 

[oo8 7] m&tt&mm*. /xjv i 6coii@D{c^ 
*5ckt>\ jtR«*£6fr*aa#fli<Df#e>tts z = 6 d io 

[0 0 8 8] ttft$tl^tt?tllll«ftU:Bmaii?et3 
5 -eti[^FtUKatK£ (StreakS) fc <fc •? iHfM 

[0 0 8 9] -7? > L D V ( K\y 7^ U- If flBSHO 
fcifcifflgfc z = 6 DfcfcttSfc^ttfcgfifciiiifcfc . 
v»T»JfflKff'3feo cn£©fliB!J«\ «MWi)We*T»* 

[0 0 9 0] ai.2fCLDVT?iiS«n/'ciB}!SaP*ff 

R e = 1 1 5 5) *^fo Hi 2£ 20 
(4, ¥£)j§fi (LDV (mean) ) kgftjlg (LD 
V (rms) ) (Dftid, e fcijVf o £fcx 01 3£S 
?*«Kffl*ffiK: J: * / X;l>PfrSrtS£6l Z = 6 Dttifi 

[0 0 9 1 ] tfcHW*H«T*tt; taa©tfrHU*iI*ffi 
(S t r e a kr£) fc<fc oT3I®£gt±iL fc 0 0 1 4*3 

*tfHi sic ^n^ r nz = 6Dtc*3it^f4^iR] I P^i3i 
.hk*oMeiraflE»#i Jtoco«i«*aau asua 

a &ttfW©4>4**fli$ fe*k B*liB73lR]fc:£*©iIi« 
*H^T'{* 1 6 0 0 0&©H{ifr5> 
IM£gUH t/Co H 1 4 fcitf H 1 5 ISLTfitWSfoft 

mts x f , g it^rmmm^ 1000 40 

T*oft, LDVfC^SfriJi'jT'f*, fttft«3 0 0 0fr5> 

[0092] eii 5tew?n**fc?isisftisa»^g f 

i *J£8t 3 t, StreakSliLDV i:fS|qJ*<-& 
*(C|gLT(±, «at(f'C> (r = 0) fcfeV^T, LDVt 

*9t>**^*#JWB£ftfc. 01 4 t^^nrcWT^iB] 

SCfctfBt&Sftfc. 50 



M2 0 0 3-8 4 0 0 5 
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[0 0 9 3] %W)M%.ft : fc*Pffl.t%t, Streak 
SttLDVfclI6jtf--&"f **>©©* LDVcD^ijjijg 
i £J±^S t r e a k£©^»jjig#*Pgti*£& 
«fOBMWIBIfr5«WU Strea 
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